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Prolonged viewing of high contrast gratings alters perceived stimulus contrast, and 
produces characteristic changes in the contrast response functions of neurons in 
the primary visual cortex (V1). This is referred to as contrast adaptation. Although 
contrast adaptation has been well-studied, its underlying neural mechanisms are not 
well-understood. Therefore, we investigated contrast adaptation in mouse V1 with the goal 
of establishing a quantitative description of this phenomenon in a genetically manipulate 
animal model. One interesting aspect of contrast adaptation that has been observed 
both perceptually and in single unit studies is its specificity for the spatial and temporal 
characteristics of the stimulus. Therefore, in the present work we determined if the 
magnitude of contrast adaptation in mouse V1 neurons was dependent on the spatial 
frequency and temporal frequency of the adapting grating. We used protocols that were 
readily comparable with previous studies in cats and primates, and also a novel contrast 
ramp stimulus that characterized the spatial and temporal specificity of contrast adaptation 
simultaneously. Similar to previous work in higher mammals, we found that contrast 
adaptation was strongest when the spatial frequency and temporal frequency of the 
adapting grating matched the test stimulus. This suggests similar mechanisms underlying 
contrast adaptation across animal models and indicates that the rapidly advancing genetic 
tools available in mice could be used to provide insights into this phenomenon. 



Keywords: adaptation, mouse vision, primary visual cortex, sinusoidal gratings, pattern-specificity, 
electrophysiology, context 



INTRODUCTION 

Our perception of the world around us, and the neural activity 
underlying this experience, is strongly dependent on the recent 
stimulus history. In the visual system, there is evidence for a 
number of self-calibration mechanisms that rapidly adapt visual 
processing according to the prevailing attributes of the stim- 
ulus being viewed (Carandini, 2000). Contrast adaptation has 
been used extensively to study this form of short-term plasticity. 
In psychophysical studies, prolonged viewing of a high-contrast 
pattern can produce a perceived fading of the adapting stim- 
ulus and reduce sensitivity to low contrasts, but it can also 
improve sensitivity and discrimination around the adapting con- 
trast (Blakemore and Campbell, 1969; Greenlee and Heitger, 
1988; Foley and Chen, 1997; Abbonizio et al., 2002). Primary 
visual cortex (VI) neurons have sigmoidal contrast response 
functions when spike rate is plotted as a function of stimulus con- 
trast, and contrast adaptation has been shown to shift the most 
sensitive part of the curve toward the adapting contrast (Movshon 
and Lennie, 1979; Ohzawa et al, 1982, 1985; Sclar et al, 1989; 
Bonds, 1991; Ibbotson, 2005). A case has also been made that 
contrast adaptation (and similar processes) must be incorporated 
into models of VI to better predict the responses of real neu- 
rons to natural stimuli (Carandini et al., 2005). Thus, there is 
converging evidence that contrast adaptation is a fundamental 
process that the visual system uses to make moment-to-moment 
adjustments in its sensitivity to incoming input. 



Both psychophysical observations and single unit recording 
studies in VI indicate that contrast adaption is pattern-specific 
such that its magnitude can depend on the spatial frequency 
(SF), temporal frequency (TF), or orientation of the adapting 
and test stimuli (Blakemore et al, 1973; Vautin and Berkley, 
1977; Movshon and Lennie, 1979; Albrecht et al, 1984; Ohzawa 
et al., 1985; Saul and Cynader, 1989a,b; Snowden and Hammett, 
1996; Miiller et al, 1999). This pattern-specificity has been used 
to constrain possible mechanisms underlying contrast adapta- 
tion. For example, both psychophysical and VI data indicate 
that contrast adaptation is strongest when the SF of the adapt- 
ing stimulus matches the test stimulus (psychophysics: Blakemore 
and Campbell, 1969; Blakemore and Nachmias, 1971; Blakemore 
et al., 1973; Snowden and Hammett, 1996; neurophysiology: 
Movshon and Lennie, 1979; Ohzawa et al, 1985; Saul and 
Cynader, 1989a), but this SF specificity must develop in the cor- 
tex because contrast adaptation in the lateral geniculate nucleus 
(LGN) does not appear to be SF specific (Duong and Freeman, 
2007). 

Several cellular and circuit mechanisms have been proposed to 
play a role in contrast adaptation (for a review see Kohn, 2007), 
but understanding of the cause of contrast adaptation remains 
incomplete. Several useful genetic tools available in mice could 
provide another avenue to explore contrast adaptation, but base- 
line conditions must first be established in this species to make 
any genetic manipulation related to contrast coding interpretable. 
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Several recent studies of mouse VI have revealed similarities 
between mice and higher mammals, including tuning for spatial 
and temporal frequencies, selectivity for orientation and direc- 
tion, and the presence of simple and complex cells (Niell and 
Stryker, 2008; Gao et al., 2010; Van den Bergh et al, 2010). 
However, contrast adaptation in mouse VI has been reported in 
only two studies (Niell and Stryker, 2008; Stroud et al., 2012). 
Stroud et al. (2012) investigated the orientation specificity of 
contrast adaptation, but the other two aspects of the pattern- 
specificity of contrast adaptation that have been so important for 
linking electrophysiological studies in higher mammals to human 
psychophysical observations, namely specificity for SF and TF, 
remain unexplored. Therefore, we examined the spatiotemporal 
specificity of contrast adaptation in mouse VI using a top-up 
adaptation protocol that was comparable with previous stud- 
ies in cat and monkey (Movshon and Lennie, 1979; Duong and 
Freeman, 2007; Dhruv et al., 2011). We also used dynamic con- 
trast ramp stimuli of varying SF and TF to obtain rapid measures 
of contrast adaptation with a wide variety of adaptors (Crowder 
et al, 2008; Stroud et al, 2012). 

Mouse VI neurons showed robust contrast adaptation when 
the adapting grating matched the neuron's preferred stim- 
ulus, which confirms earlier findings (Stroud et al., 2012). 
Furthermore, in the top-up protocol contrast adaptation was 
diminished or absent when the SF or TF of the adaptor did not 
match the neuron's preference, indicating that mouse VI neu- 
rons show adaptation specificity similar to that observed in cats 
and primates. Adaptation observed in the contrast ramp exper- 
iments was also pattern-selective, but maximal adaptation often 
occurred at slightly higher-than-preferred SFs, indicating that the 
exact properties of the contrast adaptation observed depends on 
the nature of the testing protocol. 

MATERIALS AND METHODS 
ANESTHESIA AND SURGICAL PROCEDURES 

The experimental procedures reported herein conform to the 
guidelines established by the Canadian Council on Animal Care, 
which were approved by the University Committee on Laboratory 
Animals at Dalhousie University. Electrophysiological recordings 
were made from 25 adult male C57 BL/6 J mice weighing between 
20 and 30 g, which were purchased from Jackson Laboratories 
(Bar Harbor, Maine). In early experiments, mice (n = 15) were 
sedated with chlorprothixene (5 mg/kg ip; Sigma, St. Louis, MO), 
and then anesthetized with urethane (0.5-1.2 g/kg ip; Sigma). 
If needed, a small dose of ketamine (20 mg/kg ip; Wyeth) was 
given to accelerate descent to the surgical plane of anesthesia, 
and allow a tracheotomy to be performed quickly (see Moldestad 
et al., 2009 for details). Mice were left free-breathing through- 
out the experiment and a tube located in front of the mouse 
delivered oxygen (O.lL/min) to supplement room air. In later 
experiments, mice (n = 10) were sedated with chlorprothixene 
(5 mg/kg ip) and anesthetized with isoflurane delivered through 
a customized nose cone (2.5% during induction, 1.5% during 
surgery, and 0.4-1% during recording), which decreased prepa- 
ration time by eliminating the need for a tracheotomy. Gas 
anesthesia did not appear to affect the frequency of encounter- 
ing responsive units, and produced no significant differences in 



the tuning strength or selectivity of recorded units (assessed with 
discrimination indices, see Initial data analysis below; two-sample 
t-tests, p > 0.2 for all; c.f. Kaneko et al., 2012). For all mice, body 
temperature was maintained at 37.5°C with a heating pad, and 
their corneas were protected by frequent application of a thin 
layer of optically neutral silicone oil (30000 cSt; Sigma). The skull 
was stabilized in a stereotax, and a craniotomy (~1 mm 2 ) was 
made over the monocular retinotopic representation in primary 
visual cortex (~0.8mm anterior and 2.3 mm lateral to lambda; 
Paxinos and Franklin, 2001). Recordings were made using either 
glass micropipettes (2-5 |xm tip diameter, filled with 2 M NaCl) or 
carbon-fiber in glass microelectrodes (0.6-1.5 Mf2 impedance). 
Electrode depth was controlled using a micromanipulator (FHC, 
Bowdoin, ME). Extracellular signals from individual units were 
amplified (Xcell 3+, FHC) and filtered (bandpass: 50-2000 Hz) 
before being digitized (Cambridge Electronic Design Powerl401 
with Spike2, Cambridge, England). Acquired signals were sam- 
pled at 40 kHz, and online analysis was performed on triggered 
TTL pulses with Spike2, but subsequent analysis was done offline. 

VISUAL STIMULI 

Upon isolation of a visually responsive unit, the receptive 
field (RF) was mapped using hand-driven light bars and 
spots. Quantitative testing was then performed with custom 
computer generated visual stimuli programmed in MatLab 
(MathWorks, Natick, MA) using the Psychophysics Toolbox 
extension (Brainard, 1997; Pelli, 1997), and presented on a cali- 
brated CRT monitor (LG Flatron 915FT plus 19" display, 100 Hz 
refresh, 1024 x 768 pixels, mean luminance = 30cd/m 2 ) at a 
viewing distance of 10-25 cm. All stimuli were presented in a cir- 
cular aperture surrounded by a gray field of mean luminance. 
Orientation selectivity and surround suppression were character- 
ized online using drifting square wave gratings. Spatiotemporal 
tuning was then assessed with full contrast drifting sine wave grat- 
ings with 36 combinations of SFs [0.01, 0.02, 0.04, 0.08, 0.16 and 
0.32 cycles per degree (cpd)] and TFs (0.25, 0.5, 1, 2, 4, 8 Hz). All 
spatiotemporal and adapting stimuli were presented at the opti- 
mal orientation and size for each unit, and drifted in the direction 
that elicited maximal excitation. Presentations of each combina- 
tion of SF and TF were randomized with 8-10 repeats for each 
stimulus. The presentation time of the stimulus was 1.5 s, and 
a gray of mean luminance was shown between stimuli for 0.5 s. 
Grating start-phase was staggered on each repetition to average 
out periodic firing of phase-sensitive neurons. The spatiotempo- 
ral tuning of each unit was then examined online and appropriate 
adaptors were selected for the subsequently presented contrast 
adaptation protocols. Two stimulus protocols that have previously 
been used to investigate contrast adaptation in mice, cats, and 
primates were modified to study the spatiotemporal specificity 
of contrast adaptation in mouse VI: top-up adaptation (Sclar 
et al, 1989; Duong and Freeman, 2007; Stroud et al, 2012), and 
contrast ramps (Crowder et al., 2008; Stroud et al., 2012). 

Top-up adaptation 

We chose the top-up contrast adaptation protocol because it has 
commonly been used to study the stimulus specificity of contrast 
adaptation in higher mammals (e.g., Movshon and Lennie, 1979; 
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Duong and Freeman, 2007; Dhruv et al., 2011), which facilitates 
cross-species comparisons. Sine-wave contrast is defined as: 



Michelson contrast 



(Luminance n 



Luminance m ; n ) 



(Luminance max + Luminance m ; n ) 



(1) 



where Luminance max and Luminance m i n are the maximum 
and minimum luminances, respectively. Non-adapted contrast 
response functions were obtained by recording responses to ten 
contrasts (0.04, 0.08, 0.12, 0.16, 0.24, 0.32, 0.48, 0.64, 0.82, 1) 
presented in random order for 0.5 s tests (8-12 repetitions) inter- 
leaved with 4 s of mean luminance. Adapted contrast response 
functions were collected in blocks where: (1) the adapting grat- 
ing matched the cell's spatiotemporal peak; (2) the SF of the 
adapting grating was 1-3 octaves higher or lower than the cell's 
preferred SF; and (3) the TF of the adapting grating was set to 
8 Hz. Adaptation blocks consisted of 60 s of the adapting grat- 
ing at a contrast of 0.32 followed by 0.5 s tests (aforementioned 
contrasts for 8-12 repetitions) interleaved with 4 s adaptation 
top-ups. An adapting contrast of 0.32 was chosen because our 
previous study of contrast adaptation in mouse VI (Stroud et al., 
2012) indicated that this contrast produced reliable adaptation 
while still allowing the data to be easily fit with sigmoid curves 
(see Curve Fitting below). 

Contrast ramps 

One drawback of the top-up protocol described above is that it 
takes a long time to record even a single adapted contrast response 
function (Sclar et al., 1989; Crowder et al., 2006). Therefore, when 
exploring the SF or TF specificity of contrast adaptation, only 
a few conditions can be examined for any single cell. To more 
fully assess the nature of contrast adaptation in the spatiotem- 
poral domain, we used contrast ramp stimuli. Contrast ramps 
are dynamic contrast stimuli where the contrast of the sine wave 
grating is changed linearly on each animation frame over the 
time-course of the presentation. Importantly, these ramps are 
able to measure several key markers of contrast adaptation with 
fairly short presentation times (Crowder et al, 2008; Stroud et al., 
2012). Contrast ramp stimuli were first presented at a contrast 
of 0, and contrast was increased linearly over 2 s until it reached 
1 (rising phase). The contrast of the grating was then ramped 
back down from 1 to 0 (falling phase) over the next 2 s. Thus, 
the neuron is presented with identical contrasts in the rising and 
falling phases, but the order of presentation (i.e., temporal con- 
text) is reversed. A full screen gray of mean luminance was shown 
between ramp stimuli for 2 s. In this protocol, the spatiotempo- 
ral specificity of contrast adaptation was tested by varying the SF 
and TF of the contrast ramps using the 36 combinations of SFs 
(0.01, 0.02, 0.04, 0.08, 0.16, and 0.32 cpd) and TFs (0.25, 0.5, 1, 2, 
4, 8 Hz) that were directly comparable with the spatiotemporal 
profile obtained for each neuron. Contrast ramps with differ- 
ent spatiotemporal combinations were randomized and repeated 
8-12 times for each combination. 

We were interested in determining whether the spatiotemporal 
combination that caused maximal firing also caused maximum 
hysteresis between the rising and falling portions of the contrast 
ramp. In order to test this, we used a symmetrical contrast ramp 



procedure, which maintained the same spatiotemporal parame- 
ters for both the rising and falling phase of the ramp. We also 
collected a second type of contrast ramp from a subset of neurons 
referred to as peak-tested contrast ramps that were more directly 
comparable to the top-up protocol. In the peak-tested protocol, 
the rising phase of the contrast ramp was one of the 36 combina- 
tions of SF and TF, but the falling phase was always shown at the 
neuron's preferred SF and TF (see Results), as chosen from the 
online tuning function. 

INITIAL DATA ANALYSIS 

Spike sorting was performed offline with Spike2 software, which 
first searched for and sorted spikes using a supervised template- 
matching algorithm, and then displayed candidate spikes with a 
principle components analysis for approval. Data was exported to 
MatLab and neuronal responses were represented as spike density 
functions (SDF) with 1 kHz resolution generated by convolving 
a delta function at each spike arrival time with a Gaussian win- 
dow. For each unit, we calculated the magnitude of orientation, 
size, and spatiotemporal tuning using a discrimination index (DI) 
(DeAngelis and Uka, 2003): 



DI 



( Res PMax " Res PMin) 



((Resp Max - Resp Min ) + 2VSSE/W - M)) 



(2) 



RespMax is the neuron's max response, while RespMin is the neu- 
ron's minimum response. SSE is the sum of squared error of the 
mean, N is the total number of presentations of the stimuli, and M 
is the number of different stimuli presented. In order to classify 
cells as simple or complex, we divided the first Fourier coeffi- 
cient of a neuron's response to a grating near the spatiotemporal 
peak (Pi) by the mean time-averaged response to this grating 
CF 0 ) (Movshon et al, 1978a,b; Skottun et al., 1991). Despite some 
recent controversy (Mechler and Ringach, 2002; Crowder et al., 
2007; Henry and Hawken, 2013; Hietanen et al, 2013), the Fi/F 0 
ratio has been used to quantitatively classify simple and complex 
cells in numerous studies, and an F\ /Fq ratio less than 1 indicates 
a cell is complex. 

Curve fitting 

We used the least squares method to fit contrast response func- 
tions. Sigmoid curves (Albrecht and Hamilton, 1982) were fit to 
the mean responses from top-up contrast response functions and 
SDFs produced by contrast ramps: 



R(C) 



C U + c" 



+ M 



(3) 



where R(c,) is the amplitude of the evoked response at contrast 
q, M is the spontaneous rate, n is the exponent that determines 
the steepness of the curve, R max is the maximum elevation in 
response above the spontaneous rate, and C50 is the contrast that 
generates a response elevation of half R max . Response saturation 
was evident for almost all non-adapted top-up contrast response 
functions and rising ramp responses allowing for well constrained 
fits. When fitting adapted curves where the response to maximal 
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contrast was similar to or less than the non-adapted response but 
saturation was not evident, we assigned an upper bound on the 
adapted R max of 15% above the non-adapted R max in order to 
obtain tractable fits. 

Neuronal latency 

To examine the amount of hysteresis for each contrast 
ramp, responses were latency-corrected as previously described 
(Crowder et al., 2008; Stroud et al., 2012). Briefly, for each unit 
a response threshold was established based on the 99% cut-off 
from a Poisson distribution fitted to the spontaneous firing rate. 
Each unit's response latency was calculated as the first time the 
spiking rate in the response to gratings of optimal SF and TF 
(from the spatiotemporal tuning stimulus) exceeded the afore- 
mentioned Poisson threshold and stayed above the threshold for 
the subsequent 25 ms (Price et al., 2005). For each unit, responses 
to contrast ramps were shifted back in time by the neural latency 
then split into the rising and falling phases and re-plotted using 
units of contrast on the abscissa instead of time (which resulted in 
the falling phases of contrast ramps being flipped left-to-right). 

RESULTS 

Recordings were collected from 188 visually responsive units in 
the primary visual cortex of 25 C57BL/6 J mice. We obtained con- 
trast adaptation data from 65 units using the top-up protocol, 
and 125 units using the ramp protocol (n = 90 for symmetri- 
cal contrast ramps; n = 35 for peak-tested contrast ramps). The 
stimulus preferences of units in our sample were generally consis- 
tent with previous reports. Discrimination indices for orientation 
selectivity (0.49 ±0.1; mean ± s.d.), size tuning (0.62 ±0.1), 
and spatiotemporal selectivity (0.64 ± 0.08) were similar to those 
reported by Gao et al. (2010). Peak SFs and TFs were broadly 
distributed, with preferred SFs ranging from 0.01 to 0.18 cpd 
(mean = 0.03 cpd) and preferred TFs ranging from 0.25 to 
8 Hz (mean = 1.77 Hz). Figure 1A shows the grid-like array of 
responses used to measure the spatiotemporal tuning of a sample 
neuron, and Figure IB shows how these responses can be summa- 
rized as a contour plot to indicate the combination of SF and TF 
that produced the maximal response. Our range of peak SFs and 
TFs were similar to recent electrophysiological studies of mouse 
visual cortex (Niell and Stryker, 2008; Gao et al, 2010; LeDue 
et al, 2012), and within the ranges shown by recent multi-photon 
calcium imaging studies (Andermann et al., 2011; Marshel et al., 
2011). Finally, 157 units were classified as complex (Pi /Fq ratio < 
1) and 35 units were classified as simple (F\/F$ ratio > 1). Since 
simple and complex cells showed similar trends for all measures 
of contrast adaptation, they were pooled into a single group. 

TOP-UP CONTRAST ADAPTATION 

Robust contrast adaptation following prolonged exposure to an 
adaptor of the preferred SF and TF has been shown previously in 
mouse VI (Stroud et al., 2012). However, to our knowledge the 
spatial and temporal frequency specificity of contrast adaptation 
in mouse VI have not been explored. Therefore, we compared 
the magnitude of contrast adaptation induced by an adaptor with 
preferred SF and TF with that induced by an adaptor with non- 
preferred SF or TF. We chose a non-preferred adapting TF of 8 Hz 



because high TFs rarely elicited strong responses. This high TF 
adaptor also permitted comparisons with primate work, which 
has shown that high TFs can reliably induce contrast adapta- 
tion in VI without strongly driving the recorded neurons (Dhruv 
et al., 2011), presumably by inducing adaptation in the LGN 
(Solomon et al., 2004). We selected non-preferred SFs 1-3 octaves 
higher or lower than the peak SF depending on the breadth and 
location of the recorded unit's spatiotemporal tuning. Care was 
taken to ensure non-preferred SFs elicited weak responses from 
the recorded unit but also were within the range of peak SFs of 
our sample population and below the mean SF cutoff reported in 
previous studies of LGN and VI (Grubb and Thompson, 2003; 
Gao et al., 2010). Higher adapting SFs were selected more often 
than lower ones since SFs lower than 0.01 cpd can begin to appear 
as global changes in luminance within the stimulus aperture. 

Figures 1C-F shows the SF and TF specificity of contrast adap- 
tation for four example neurons. Contrast response functions 
are shown for non-adapted (black squares), preferred adapted 
(red circles), non-preferred SF adapted (green triangles), and 
non-preferred TF adapted (blue stars) conditions. For the cell 
in Figure IF, contrast response functions from two different 
non-preferred SFs (low SF = pink diamonds; high SF = green 
triangles) are shown. The spatiotemporal tuning of each unit is 
shown inset with the SF and TF values of the adapting stim- 
uli indicated with matching symbols. In each case the preferred 
adaptor induced the most contrast adaptation. Non-preferred 
adaptors either induced virtually no adaptation (Figures 1C,D), 
or less adaptation than the preferred stimulus (Figures 1E,F). 

Sigmoid fits to each contrast response function are shown 
as thin lines in Figure 1, and we used the C50 and Rmax 
parameters extracted from these fits to quantitatively analyze 
changes in contrast response functions following top-up adap- 
tation. For each adaptation condition we measured the change 
from the non-adapted curve as a difference-over-sum calcula- 
tion (parameter sm ft = [adapted -non-adapted]/ [adapted + non- 
adapted]), and plotted this metric as population histograms in 
Figure 2. For Figures 2A,C,E positive values of C5o- s hift indicate a 
rightward shift in the adapted contrast response function. Nearly 
all cells showed a rightward shift following preferred adapta- 
tion (Figure 2A, mean cso-shift = 0.26), but the population was 
centered closer to zero for both adaptors with non-preferred 
SF (Figure 2C, mean c 50 -shift = 0.09) and TF (Figure 2E, mean 
cso-shift = 0.05). A One-Way ANOVA followed by a Tukey- 
Kramer post-hoc indicated that the preferred adaptation pro- 
duced significantly larger values of cso-shiftthan the other two 
adaptation conditions [F(2, 169) = 27.54, p < 0.001], while non- 
preferred SF and TF Cso-shift did not differ. For Figures 2B,D,F 
negative values of R max -shift indicate a decrease in firing to max- 
imal contrast following adaptation. Most cells showed a mod- 
est decrease in R max following preferred adaptation (Figure 2B, 
mean R m ax-shift = —0.17), but the population was centered near 
zero for both adaptors with non-preferred SF (Figure 2D, mean 
Rmax-shift = 0.02) and TF (Figure 2F, mean R max - sm ft = 0.05). A 
One-Way ANOVA followed by a Tukey- Kramer post-hoc showed 
similar results to the C50 data, with preferred adaptation pro- 
ducing significantly more negative values of R max - S hift than the 
other 2 adaptation conditions [Fp, 169) = 7.72, p < 0.001], while 
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FIGURE 1 | Spatiotemporal specificity of top-up contrast adaptation. 

The spatiotemporal selectivity of a sample neuron is presented in (A), with 
a grid of SDFs showing the neuron's response to gratings with different 
combinations of SF (columns) and TF (rows). A scale bar depicting time 
vs. impulses per second (ips) is shown in the bottom right SDF (0.32 cpd 
and 0.25 Hz). Mean responses from SDFs are summarized as a grayscale 
contour plot in (B) (maximal firing shown in white), with SF on the 
abscissa and TF ordinate. (C-F) show contrast response functions from 
four sample neurons with contrast on the abscissa and mean response 



rate on the ordinate. Non-adapted, preferred adapted, non-preferred SF 
adapted, and non-preferred TF adapted responses are shown as black 
squares, red circles, green triangles, and blue pentagrams, respectively. (F) 
shows responses following an additional non-preferred SF adaptor as pink 
diamonds. Error bars represent SEM, and thin lines represent best fits to a 
sigmoid function (see Materials and Methods). Spatiotemporal tuning of 
each sample neuron is shown in the corresponding inset with adapting 
stimulus SF and TF indicated with symbols matched to contrast response 
functions. 



non-preferred SF and TF R max -shift did not differ. Another way of 
quantifying the spatiotemporal specificity of contrast adaptation 
is simply to rank order the adapted curves for each cell. Preferred 
adaptation C50 values were larger than C50 values measured fol- 
lowing non-preferred SF adaptation for 90% of cells, and non- 
preferred TF adaptation for 92% of cells. Preferred adaptation 
Rmax values were smaller than R max values measured following 
non-preferred SF adaptation for 75% of cells, and non-preferred 
TF adaptation for 70% of cells. 

CONTRAST RAMP ADAPTATION 

The top-up adaptation data above clearly demonstrates that the 
magnitude of contrast adaptation in mouse VI depends on the 
adapting SF and TF, however, as noted in the Methods section 
only a few adaptation conditions can be studied for any sin- 
gle cell due to the time constraints imposed by this protocol. 



Therefore, we used symmetrical contrast ramp stimuli to more 
extensively map the spatiotemporal selectivity of contrast adap- 
tation. Figure 3A shows the response of a representative neuron 
to a contrast ramp of optimal SF and TF. Even though the ris- 
ing and falling phases of the ramp stimulus are symmetrical, the 
spiking response shows clear hysteresis. If this spiking response is 
latency-corrected and re-plotted with contrast on the abscissa (see 
Materials and Methods), the difference between the responses to 
the rising (red) and falling (blue) phases of the contrast ramp is 
accentuated further (Figure 3B). As in previous studies (Crowder 
et al, 2008; Stroud et al., 2012), the SDFs were fit to sigmoid 
curves (thin lines). The most useful parameter extracted from 
the sigmoid fits was C50, since it captured the rightward shift 
in the contrast response function by comparing semi-saturation 
contrasts of the rising (upward pointing arrowhead) and falling 
phases (downward pointing arrowhead) of the contrast ramp. For 
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FIGURE 2 | Population data from top-up adaptation. (A,C,E) plot c 5 o-shift 
population histograms following preferred adaptation, adaptation with 
non-preferred SFs, and adaptation with non-preferred TFs, respectively. 
(B,D,F) plot Fimax-shift population histograms following preferred adaptation, 
adaptation with non-preferred SFs, and non-preferred TFs, respectively. 
Arrowheads represent population means, and asterisks denote significant 
differences (see Results). 
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FIGURE 3 | Contrast ramps of optimal SF and TF. (A) shows the 
response of a representative neuron to a contrast ramp of optimal SF and 
TF along with a schematic demarcating the rising and falling phases of the 
contrast ramp with dashed lines. In (B), the neuron's responses to the 
rising (blue lines) and falling phases (red lines) of the contrast ramp are 
compared by folding the latency-corrected SDF back on itself so that 
contrast is on the abscissa and spikes/s is on the ordinate (see Materials 
and Methods). SDFs were fit with sigmoid functions (thin lines), and 
upward and downward pointing arrowheads represent c 5 o values obtained 
from fits to the rising and falling phases, respectively. Population data 
comparing c 5 o values obtained from fits to the rising (abscissa) and falling 
phase responses (ordinate) is shown in (C). 



our sample (n = 90), C50 values from the rising phase were almost 
always smaller than C50S from the falling phase (Figure 3C), and 
this difference was significant (p < 0.01, paired f-test). This repli- 
cates earlier findings in cats (Crowder et al., 2008), and mice 
(Stroud etal.,2012). 

To map the spatiotemporal specificity of contrast adaptation 
we measured the hysteresis of ramp responses when the SF and 
TF of the ramp grating were varied (for easy comparison to the 
spatiotemporal tuning also obtained for each neuron we used 
the same 36 combinations of SF and TF). This stimulus protocol 
examined the spatiotemporal specificity of contrast adaptation 
from a slightly different perspective than the top-up protocol. 
The top-up protocol measured whether the magnitude of contrast 
adaptation was affected if the adapting grating did not match the 
test grating, which emphasized the importance of the adapting 
stimulus. Symmetrical contrast ramps measured the combina- 
tion of SF and TF that produced the most hysteresis, which 
emphasized the importance of the cell's own preferred stimu- 
lus in determining the strength and specificity of the adaptation 
effect (Saul and Cynader, 1989a). Figure 4A shows a grid of SDF 
ramp responses from a sample cell, each with the same format 
as Figure 3B. This neuron had strong ramp responses with sub- 
stantial hysteresis around 0.02-0.04 cpd and 1-2 Hz. Responses to 
lower TFs (~0.25 Hz) showed little hysteresis despite monotonic 
increases in firing with contrast, and the entire ramp response 



flattened out at the highest SFs and TFs. The former effect was 
observed in 81/90 neurons, indicating that diminished adapta- 
tion was not solely due to lack of responding. We wanted to 
summarize the pattern of hysteresis for each neuron as a con- 
tour plot, but C50 values taken from sigmoid fits to SDFs were 
unreliable for spatiotemporal combinations away from the peak, 
so we measured adaptation by calculating the mean difference 
between the responses to rising and falling phases of the con- 
trast ramps. During adaptation, the semi-saturation contrast of 
the falling phase ramp response shifts to higher values, causing the 
falling ramp response to be lower than the rising phase response 
at most contrasts. This method of analysis has been shown by 
Stroud et al. (2012) to capture the major features of contrast 
adaptation without relying on fitting the ramp SDFs to sigmoid 
functions. Figure 4B shows the contour plot for this neuron sum- 
marizing the magnitude of hysteresis evoked by each combination 
of SF and TF. The first feature to note is the clear peak around 
0.04 cpd and 1-2 Hz. Likewise, 85 out of 90 units produced con- 
tour plots with an easily identifiable single peak that was at least 
four times higher than the level of hysteresis produced by the least 
effective ramp. This supports our earlier finding of spatiotempo- 
ral specificity of contrast adaptation using a different method. 
The second feature of the hysteresis contour plot that we were 
interested in was whether the combination of SF and TF that 
produced maximum hysteresis for contrast ramps matched the 
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FIGURE 4 | Spatiotemporal specificity of contrast ramp adaptation. 

The grid of SDFs in (A) show the hysteresis induced in a sample neuron by 
symmetrical contrast ramps of varying TFs (rows) and SFs (columns). Each 
SDF follows a similar format to Figure 3B, with responses to the rising 
phase of the contrast ramp shown in blue and responses to the falling 
phase shown in red. A scale bar depicting time vs. impulses per second 
(ips) is shown in the lower right (0.32 cpd and 0.5 Hz). The spatiotemporal 
pattern of hysteresis for the sample neuron is represented as a blue-tinted 
contour plot in (B), with larger mean differences between the rising and 
falling phase responses shown as more desaturated hues (see Results). For 
comparison, the spatiotemporal tuning of the sample neuron is shown in 
(C) as a grayscale contour plot, and the correlation between the two 
contour plots is indicated (double-headed arrow). For both contour plots SF 
is on the abscissa and TF in on the ordinate. 



neuron's peak in the spatiotemporal domain tested with regular 
grating blocks (Figure 4C). For this neuron, the two contour plots 
look similar (R = 0.78 from a 2D correlation analysis), but the 
gratings that produced maximum hysteresis had a slightly higher 
SF than the gratings that produced maximum firing. Figure 5 
shows two more example cells, one where the spatiotemporal 
locations of maximum firing and maximum hysteresis match 
quite closely (Figures 5A,B; R = 0.87), and another where spa- 
tiotemporal location of maximum hysteresis is at a higher SF and 
lower TF (Figures 5C,D; R = 0.51). Figure 5E plots the differ- 
ence in peak locations from the two types of contour plots in 
the spatiotemporal domain for each cell as "hatpins" (« = 85), 
with the empty dots indicating the location of maximum hys- 
teresis. No pattern was apparent, indicating that there was not 
one specific combination of SF and TF that universally induced 
maximal hysteresis across cells. Figure 5F normalized the data 
from Figure 5E by calculating the octave difference in SF and 



TF between each pair of peaks to show the location of maxi- 
mal contrast hysteresis (empty dots) relative to each cell's peak in 
spatiotemporal tuning (all normalized to 0). Although the differ- 
ences were small (45% of cells had peaks within 1 octave of each 
other, and the median .R-value from 2D correlations was 0.71), at 
the population level the gratings that produced maximum hys- 
teresis tended to have slightly higher SFs (mean: 0.47 octaves) 
and lower TFs (mean = —0.22 octaves) than the gratings that 
produced maximum firing. A 2-way repeated measures ANOVA 
(with peaks from spatiotemporal tuning vs. hysteresis contour 
plots and SF vs. TF as factors) showed a significant main effect 
of peak type [F(i_ 84) = 3.95, p < 0.05], indicating that the spa- 
tiotemporal location of peak hysteresis and peak firing tended to 
be different. Furthermore, a significant interaction between fac- 
tors indicated that the difference between peaks was larger for SF 
than for TF [F (h 84) = 15.75, p < 0.001]. 

Considering that we consistently observed the strongest adap- 
tation when using the preferred SF and TF in the top-up protocol 
we were surprised by the results of the ramp protocol. However, as 
noted above there was one key difference between adapting proce- 
dures: in the top-up protocol the adapting grating varied but the 
test gratings were always set at the preferred SF and TF, whereas 
in the symmetrical ramp protocol the SF and TF of the grating 
remained constant throughout the rising and falling phases of 
the contrast ramp. To determine whether the different adaptation 
effects observed between the top-up and contrast ramp proto- 
cols were due to switching the SF/TF between adapting and test 
stimuli or the dynamic nature of the contrast ramp stimuli we 
presented a subset of cells with a modified ramp protocol referred 
to as peak-tested ramps. For these peak- tested ramps, the rising 
phase could have any one of the 36 combinations of spatiotem- 
poral frequencies (a proxy for the adapting gratings in the top-up 
protocol), but the falling phase was always shown at the neuron's 
preferred spatiotemporal frequency (a proxy for the test gratings). 
For these stimuli, we compared the response to the falling phase at 
the spatiotemporal peak with the falling phase responses at every 
other spatiotemporal combination since these stimuli were iden- 
tical with only the preceding rising phase differing (Figure 6A). 
We expected the difference to be large if no contrast adaptation 
occurred, or small if contrast adaptation did occur. Responses 
from a representative neuron are shown in Figures 6A-C. We 
again represented the spatiotemporal specificity of adaptation 
for each neuron as a contour plot (Figure 6B), compared the 
adaptation contour plot to each neuron's spatiotemporal profile 
(Figure 6C), and calculated the octave difference in SF and TF 
between peaks for the population (Figure 6D). For peak-tested 
ramps, there were clear peaks in the hysteresis contour plots of 
every cell (e.g., Figure 6B). Importantly, the contrast hystere- 
sis and spatiotemporal profile contour plots were much more 
similar using this protocol. For the sample neuron shown in 
Figures 6A-C the 2D correlation between contour plots was 0.96, 
and the population median was 0.86. Furthermore, the 2D cor- 
relations between contour plots for the peak-tested protocol were 
significantly higher than for the symmetrical ramp protocol (p < 
0.0001; f-test). Mean octave differences in peak location between 
contour plots were —0.08 and —0.06 for SF and TF, respectively 
(Figure 6D). A Two-Way repeated measures ANOVA (with peaks 
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FIGURE 5 | Population data from symmetrical contrast ramp 
adaptation. The spatiotemporal specificity of adaptation induced with 
contrast ramps (left column) and spatiotemporal tuning (right column) are 
compared for two additional sample neurons. As in Figure 4, contrast ramp 
hysteresis is shown as blue-tinted contour plots, and spatiotemporal tuning 
is shown as grayscale contour plots. The contour plots from the neuron in 
the top row (A,B) match quite closely, indicating that similar grating 
parameters produced maximum firing and maximum contrast ramp 
hysteresis. The contour plots from the second neuron (middle row; C,D) 
match less closely, with the spatiotemporal location of maximum hysteresis 
occurring at a higher SF and lower TF than the peak in spatiotemporal 
tuning. The correlations between the contrast ramp and spatiotemporal 
tuning contour plots are indicated for each neuron (double-headed arrows). 
(E) shows population data comparing the spatiotemporal location of 
maximal hysteresis (empty dots) with the locations of peak responding 
from spatiotemporal tuning (lines). For (A-E), SF is on the abscissa and TF 
in on the ordinate. (F) normalized the data from (E) by calculating octave 
differences in SF (abscissa) and TF (ordinate) to show the location of 
maximal contrast hysteresis (empty dots) relative to each cell's 
spatiotemporal tuning (all normalized to 0). Population mean is shown as a 
solid red circle. 



from spatiotemporal tuning vs. hysteresis contour plots and SF 
vs. TF as factors) indicated that neither the main effect of peak 
type [F(i t 34) = 0.31, p > 0.57], nor the interaction between fac- 
tors were significant [F^ 34) = 0.01, p > 0.92]. Thus, switching 
the SF/TF between adapting and test stimuli appear to be the 
important difference between the top-up and symmetrical ramp 
protocols because when the contrast ramp stimulus was altered to 
more closely resemble the top-up protocol the adaptation effects 
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FIGURE 6 | Peak-tested contrast ramp adaptation. The grid of SDFs in 
(A) show the responses of a sample neuron to peak-tested ramps where 
the TFs (rows) and SFs (columns) of the rising phase of the ramp were 
varied, but the falling phase was always shown at the neuron's peak SF and 
TF (red lines). The transition between non-preferred and preferred gratings 
is especially apparent at high SFs. A scale bar depicting time vs. impulses 
per second (ips) is shown in the top left (0.01 cpd and 8 Hz). The 
spatiotemporal pattern of adaptation for the sample neuron is represented 
as a red-tinted contour plot in (B), with smaller mean differences between 
the preferred falling phase and other falling phase responses shown as 
more desaturated hues (see Results). For comparison, the spatiotemporal 
tuning of the sample neuron is shown in (C) as a grayscale contour plot, 
and the correlation between the two contour plots is indicated 
(double-headed arrow). For both contour plots SF is on the abscissa and TF 
in on the ordinate. (D) shows the octave differences in SF (abscissa) and TF 
(ordinate) between the locations of maximal contrast adaptation (empty 
dots) relative to each cell's spatiotemporal tuning (all normalized to 0). 
Population mean is shown as a solid pink circle. 



also matched the top-up results more closely. Overall, each of 
the three adaptation protocols along with their differing methods 
of analysis demonstrated the spatiotemporal specificity of con- 
trast adaptation in mouse VI, even though differences between 
protocols produced some subtle variations in the nature of the 
adaptation. 

DISCUSSION 

This study demonstrated that contrast adaptation in mouse VI 
is specific in the spatiotemporal domain. The magnitude of con- 
trast adaptation observed in single units was found to depend on 
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both the SF and TF of the adapting grating, and the nature of the 
adaptation effect could also be affected by the SF and TF of the 
test stimuli. The properties of contrast adaptation we observed 
were broadly similar to single unit studies in higher mammals 
(monkeys: Sclar et al., 1989; Dhruv et al, 2011; cats: Movshon 
and Lennie, 1979; Ohzawa et al., 1982, 1985; Saul and Cynader, 
1989a,b; Bonds, 1991) and psychophysical data (e.g., Blakemore 
and Campbell, 1969). This suggests that contrast adaptation can 
be thought of as a general feature of the mammalian geniculo- 
striate pathway along with other classical response properties 
(e.g., Niell and Stryker, 2008; Gao et al, 2010; Van den Bergh 
et al., 2010). Despite marked differences between animal mod- 
els (frontal eyes vs. lateral eyes, nocturnal vs. diurnal, acuity that 
varies over several orders of magnitude), adaptation in mouse 
visual cortex appears to follow similar rules and is of similar com- 
plexity to higher mammals. We believe that these findings uphold 
the viability of the mouse model for studying vision, and support 
the validity of a multi-species approach for investigating cortical 
visual processing. 

COMPARISON WITH PREVIOUS STUDIES 

To our knowledge only two previous studies have investigated 
contrast adaptation in mouse VI (Niell and Stryker, 2008; Stroud 
et al., 2012). Stroud et al. (2012) were able to make direct com- 
parisons between adaptation in mouse and cat VI neurons, and 
reported that most key features of contrast adaptation were sim- 
ilar between species. When adapted and tested with an optimal 
grating, adaptation shifted contrast response functions down and 
to the right. Moreover, contrast ramps produced relatively robust 
contrast adaptation given their brief presentation times. The cur- 
rent study is in agreement with these previous findings, so the 
Discussion will focus on the spatiotemporal specificity of contrast 
adaptation. 

Several previous papers have used some version of the top-up 
protocol to investigate either the SF or TF dependence of con- 
trast adaptation (e.g., Movshon and Lennie, 1979; Duong and 
Freeman, 2007; Dhruv et al., 201 1), and are therefore readily com- 
parable to our own top-up data. Within our top-up protocol, the 
test stimuli were always at the optimal SF and TF for each neuron, 
which is most similar to the stimuli used by Dhruv et al. (201 1 ) in 
their study of TF and orientation specificity of contrast adaptation 
in macaque VI. For these two studies, one central question was 
whether an adapting stimulus that itself does not strongly drive 
the recorded neuron could induce adaptation. Another study 
examining the SF specificity of contrast adaptation in cat VI used 
a complementary design where the test gratings were not at each 
neuron's optimal SF, but rather at an SF that evoked approxi- 
mately the same firing rate as the adapting grating (Movshon 
and Lennie, 1979). Regardless of these design differences, the gen- 
eral findings of these studies indicate that contrast adaptation is 
most robust when the parameters of the adapting grating are sim- 
ilar to the test grating. The current study extends this finding 
into a genetically tractable animal model where there exists an 
expanded toolbox to investigate the mechanisms underlying this 
specificity. 

Comparing orientation specificity and spatiotemporal speci- 
ficity of contrast adaptation in mouse VI is also worthwhile. Most 



neurons in VI adapt to any orientation, even ones that elicit low 
firing rates (Stroud et al, 2012), which contrasts with our current 
results in the spatiotemporal domain. This pattern of results could 
be produced if cortical adaptation mechanisms pooled over ori- 
entation (or the sharpness of tuning was diluted by non-oriented 
cells), but were at least somewhat selective in the spatiotemporal 
domain (Andermann et al., 2011; LeDue et al., 2012). 

It has been shown that adapting gratings with a high TF 
can induce modest but reliable contrast adaptation in macaque 
VI without strongly driving the recorded neuron (Dhruv et al., 
201 1), presumably by inducing adaptation in magnocellular cells 
in the LGN (Solomon et al, 2004). Therefore, we were somewhat 
surprised to observe only occasional adaptation to higher TFs in 
our data set. In their study of macaque VI, Dhruv et al. (2011) 
used an adaptor with a TF of 30-50 Hz, which was 2-3 times 
higher than the peakTF of LGN neurons (10-16 Hz: Derrington 
and Lennie, 1984; Hawken et al., 1996). For our top-up protocol, 
the high TF adaptor (8 Hz) was also approximately double the 
peak TF of mouse LGN neurons (3.8 Hz: Grubb and Thompson, 
2003), yet we observed little consistent adaptation to this stimu- 
lus. We had also predicted that the peaks in the contrast hysteresis 
contour plots obtained with our contrast ramp stimuli may be 
skewed toward higher TFs, but this was not the case. We are 
unsure what underlies this apparent species difference, but as out- 
lined in a model of multiple sources of adaptation used by Dhruv 
et al. (2011), it suggests that less adaptation is occurring in (or 
being inherited from) the LGN in mice. This observation, in con- 
junction with the finding that contrast adaptation in cat LGN 
does not show SF specificity (Duong and Freeman, 2007), pro- 
vide two good reasons for future work to investigate adaptation 
in mouse LGN. 

Finally, the contrast ramp stimuli we used in this study were 
quite unique and therefore less comparable to previous papers 
(although the staircase-like stimuli used by Bonds (1991) also 
measured hysteresis when contrasts were presented in an ordered 
manner). However, this data is relevant to a longstanding issue 
in the contrast adaptation literature. Vautin and Berkley (1977) 
were the first to discuss how the adaptation measured in a 
recorded neuron could arise from processes occurring within 
the cell itself (intrinsic) or be inherited from other neurons in 
the circuit/network (extrinsic). Saul and Cynader (1989a) sug- 
gested that intrinsic mechanisms may be more narrowly tuned 
since they depend on the cell's own tuning, while the aforemen- 
tioned model by Dhruv et al. (201 1) specified that some extrinsic 
sources of adaptation should be broadly tuned for certain stimu- 
lus attributes like orientation. The main strength of the contrast 
ramp stimulus is that it can measure features of contrast adap- 
tation on a relatively short time-scale, and this allows for a large 
stimulus-space to be explored in a reasonable amount of time. 
This seems ideal for exploring the putative differences in tun- 
ing between intrinsic and extrinsic sources of adaptation. In the 
current study, both the symmetrical and peak-tested contrast 
ramp protocols support the spatiotemporal specificity of contrast 
adaptation initially described using the top-up protocol despite 
the fact ramp stimuli measured adaptation on a different time 
scale and used different metrics. It would be interesting to obtain 
comparative data from higher mammals for these stimuli. 
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USING MOUSE MODELS TO STUDY CONTRAST ADAPTATION 

Electrophysiological studies in various animal models additively 
suggest that contrast adaptation is initiated at pre-cortical stages 
(retina: reviewed in Demb, 2008; LGN: Sanches-Vives et al., 
2000a; Solomon et al., 2004; Duong and Freeman, 2007), and 
then refined and strengthened in the cortex (Ohzawa et al., 1985; 
Carandini, 2000; Dhruv et al., 2011). Furthermore, hyperpo- 
larization of the membrane potential has been associated with 
contrast adaptation (Carandini and Ferster, 1997; Sanches-Vives 
et al, 2000a,b). However, questions about the specific mecha- 
nisms involved, their relative contributions, and the stage each 
one is implemented remain unanswered. It seems that some of 
the investigative tools currently most readily applied in the mouse 
could provide insights into these cellular and circuit mechanisms. 
The same biochemical and genetic flexibility that has allowed the 
use of optogenetic modulation to attribute particular functions 
to genetically defined inhibitory neurons within mouse VI (e.g., 
Adesnik et al., 2012; Atallah et al., 2012), or allowed genetically 
encoded calcium-indicator proteins to explore the response prop- 
erties of hundreds of visually responsive neurons simultaneously 
(e.g., Andermann et al., 2011), could also be used to explore the 
mechanisms underlying contrast adaptation. Moreover, if specific 
mechanisms are isolated they could be knocked-out or modulated 



in real-time to probe the perceptual relevance of contrast adap- 
tation using psychophysical tasks developed for the mouse (e.g., 
Busse et al., 2011). The possibility of causally linking neural 
processing to contrast perception is especially intriguing consid- 
ering psychophysical studies of the performance enhancement 
conferred by contrast adaptation have been somewhat equivocal 
(Barlow et al., 1976; Maattanen and Koenderink, 1991; Abbonizio 
et al, 2002; Kohn, 2007). We are hopeful that insights gleaned 
from the mouse model will be relevant to higher mammals 
because the properties of cortical contrast adaptation that have 
already been explored appear quite similar between species. 
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